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Abstract. We compare two predictions regarding the microscopic fluctuations of 
a system that is driven away from equilibrium: one due to Crooks [J. Stat. Phys. 
90, 1481 (1998)] which has gained recent attention in the context of nonequilibrium 
work and fluctuation theorems, and an earlier, analogous result obtained by Bochkov 
and Kuzovlev [Zh. Eksp. Teor. Fiz. 72(1), 238247 (1977)]. Both results quantify 
irreversible behavior by comparing probabilities of observing particular microscopic 
trajectories during thermodynamic processes related by time-reversal, and both are 
expressed in terms of the work performed when driving the system away from 
equilibrium. By deriving these two predictions within a single, Hamiltonian framework, 
we clarify the precise relationship between them, and discuss how the different 
definitions of work used by the two sets of authors gives rise to different physical 
interpretations. We then obtain a extended fluctuation relation that contains both the 
Crooks and the Bochkov-Kuzovlev results as special cases. 
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1. Introduction 

Recent interest in the nonequilibrium thermodynamics of small systems [Ij has brought 
attention to a set of papers written by Bochkov and Kuzovlev in the late 1970's and 
early 1980's [21 El HI [5]. In view of evident similarities between results derived in these 
papers, and predictions obtained over the past decade or so by other authors (see in 
particular [6l [TJ [H [9l [lO]), it is desirable to clarify the precise relationships between 
between these two sets of results. In an earlier paper [TT], one of us has compared the 
nonequilibrium work relation of [6], (e~^'^) = e~^^^, with an analogous equality found 
in [2], [3l m |5] , (e"^^") = 1. In the present paper we turn our attention to fluctuation 
relations for systems driven away from equilibrium. Specifically, we compare a result 
due to Crooks [8| (([1]) below) with a similar result obtained by Bochkov and Kuzovlev 

Briefly, the setting is a system driven away from equilibrium by an externally 
applied, time-dependent force, Xf. We consider two different protocols (or schedules) 
for applying this force, Xf and X^, related by time- reversal ([5]), where the superscripts 
denote "forward" and "reverse" . The microscopic evolution of the system in either case 
is described by a phase space trajectory, 7. Under appropriate assumptions regarding 
this evolution. Crooks has shown that the probability to observe a given trajectory 7^ 
during the forward process, and the probability to observe its time-reversed counterpart 
7"^ during the reverse process, are related by [8] 

_ f3(W-AF) /-,^ 

Here, W denotes the work performed on the system, and AF is a free energy difference 
between the equilibrium ensembles from which initial conditions are sampled, as 
described in detail below. From this result others are easily derived, including a 
fluctuation theorem for the ratio of forward and reverse probability distributions of work 
values, Pf{W) and Pr(W) [9], which has recently been confirmed experimentally [T2] . 

In a similar context, Bochkov and Kuzovlev have made an analogous assertion 
(Equation 2.11 of [4j or equation 7 of [3]), which in the notation of the present paper 
reads, 



^mo^ (2) 



P0[7«] 

where Wo again denotes work. 

As the above notation suggests, the two definitions of work, W and Wq, differ, 
as do the probability distributions Pf/r and Pp/j^- Moreover, AF appears in ([T]), but 
not in ([2]). To establish the relationship between these two predictions, we first derive 
them both within a Hamiltonian framework, in section O Following that, in section [3] 
we derive a more general result, containing both ([1]) and ([2]) as special cases. This 
extended relation is obtained for both Hamiltonian dynamics and Markovian stochastic 
dynamics. 



Comparison of work fluctuation relations 



3 



Independently, Seifert [IB] has concluded that both ([T]) and ([2D are special cases 
of a more general result, obtained within a theoretical framework that defines entropy 
production in out-of-equilibrium systems [H]. Our approach complements Seifert's, by 
focusing on the different physical interpretations of work appearing in and ([2]). 

2. Derivations of Earlier Results 

In this section we derive ([1]) and ([2]) within a single framework, using Hamilton's 
equations to model the microscopic evolution of the system of interest. We introduce 
and specify the notation for this framework in section 12.11 and then in section 12.21 and 
12.31 we obtain the desired results. These side-by-side derivations clarify the physical 
interpretation of the differences between ([1]) and ([2]), as highlighted in a brief summary 
at the end of this section. We note that Crooks originally obtained ([T]) by modeling the 
evolution of the system as a stochastic, Markov process. The Hamiltonian derivation of 
this result presented here is similar to the analyses of [151 IIH] • 

2.1. Setup 

Consider a classical system with N degrees of freedom, described by coordinates 
q = (gi, . . . , qn) and conjugate momenta p = {pi, . . . ,Pn), and let z = (q, p) denote a 
point in its phase space. We will be interested in the evolution of this system, in the 
presence of an externally controlled, time-dependent force Xt. In this section we model 
this evolution using Hamilton's equations, assuming a time-dependent Hamiltonian of 
the form 

H{z;Xt)=Ho{z)-Xta{z). (3) 

Here a is the coordinate conjugate to the external force X, and the "bare" Hamiltonian 
Hq denotes the energy of the system in the absence of this force. For simplicity, we 
further assume that H is time-reversal-invariant for any fixed value of the external 
force, that is, 

Hiz*;X)=Hiz;X), (4) 

where the asterisk denotes a reversal of momenta, p —p. 

We have assumed that X appears linearly in the definition of the Hamiltonian, 
H = Hq — Xa, as in the papers by Bochkov and Kuzovlev [21 [3l HI [5] . As mentioned by 
those authors, this assumption can be relaxed to allow for nonlinear coupling H{z;X), 
as in [SI [H m El Ho] . However, for clarity of exposition, we follow [21 El IH [5] and stay 
with the linear assumption throughout this paper. 

We will use the term process to denote the following sequence of events. First, the 
system is prepared in a state of thermal equilibrium at inverse temperature p. (This step 
will be discussed in more detail below.) Subsequently, from t = to t = t the system 
evolves in time under Hamilton's equations, as the external force is applied according 
to a schedule, or protocol, Xt. We will use the notation zt to specify the phase space 
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Figure 1. Schematic depiction of a forward trajectory 7^ and its conjugate twin 7^. 
The two are related by a reversal of momenta and the direction of time, i.e. = z^Zt 
for < i < T. 



coordinates of the system at a particular time t during this interval, and 7 to denote the 
entire trajectory from t = to t = r. Such a trajectory represents a single realization 
- i.e. one possible microscopic history - of the process in question. By performing this 
process repeatedly - always preparing the system in equilibrium, and always applying 
the same protocol Xt - we effectively generate independent samples 71, 72, ■ ' ' from a 
statistical ensemble of possible realizations of the process. 

We will consider two different processes, which we label forward (F) and reverse 
(R). During the forward process the force is switched from = A to X^ = B using 
the forward protocol X^ , while during the reverse process the force is switched from B 
back to y4, using the time-reversed protocol, i.e. 

= (5) 

Our assumption of time-reversal invariance then implies that solutions of Hamilton's 
equations come in conjugate pairs [Ul [IS]: for every trajectory 7-^ = {zf} that 
represents a microscopic realization of the forward process, the conjugate twin 7^ = 
{zf^} is a realization of the reverse process, where 

(6) 

as illustrated in figure [H While we restrict our analysis to Hamiltonian dynamics, 
we note that conjugate pairing of trajectories holds also for non-Hamiltonian, time- 
reversible deterministic dynamics, as exploited by Evans to derive nonequilibrium work 
relations for deterministically thermostatted systems [17] . 
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In the derivations that follow, we will make use of two different definitions of the 
work performed on the system during such a process. The first definition, which appears 
in the papers by Bochkov and Kuzovlev (e.g. equation 2.9 of [1]), is the integral of force 
versus displacement familiar from introductory mechanics texts [T8] : 

Wo = - a — dt= Xadt = Ho{zr)-Ho{zo). (7) 



JO da JO 

By contrast, recent papers (see e.g. [6], [8]) have used a different definition, 

X— dt = -y Xadt = H{zr;Xr)-H{zo;Xo) (8) 

which traces its origins to discussions of the statistical foundations of thermodynam- 
ics [ini [20l [21]. In both ([7]) and ([8]) the properties of Hamilton's equations have been 
invoked to rewrite the integral as a net change in the value of a Hamiltonian, either 
Hq or H{see [H]). At the end of this section we briefly discuss the distinction between 
these two definitions of work. 

2.2. Crooks Fluctuation Relation 

For the situation that we have just described, we now derive ([1]). In this equation, Pp['y^] 
denotes the probability distribution to observe a trajectory 7^, when performing the 
forward process. This distribution must be defined with respect to a measure on the 
space of trajectories 7^. Since Hamiltonian evolution is deterministic, it is natural to 
use the Liouville phase space measure, applied to the initial conditions of the trajectory. 
Thus, if we imagine a narrow "tube" of Hamiltonian trajectories evolving from t = 
to t = T, from a small patch of initial conditions in phase space, then we equate 
the path-space measure of this tube of trajectories with the Liouville measure of the 
patch of initial conditions: d7^ = d^;^ = d^qo d^po. Similarly, for the reverse process 
the distribution PrI'j^] is defined with respect to the measure d7^ = dzQ. Because 
Hamiltonian evolution preserves phase space volume [22], the measure of a given set of 
forward trajectories is the same as that of the conjugate set of reverse trajectories: 

d7^ = d7^. (9) 

In Crooks's formulation [8], the initial conditions for the forward process are 
sampled from the equilibrium (canonical) distribution corresponding to the initial value 
of the external force, Xq = A: 

PfIi^ = P7i4) = ^ eM-m^^; A)]. (10) 

Similarly for the reverse process, 

P^[7^] =^^^(^0"^) = ^^exp[-(3H{z^;B)]. (11) 
In these equations 

Z{X) = [ dze-^^^''^^ (12) 
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is the classical partition function associated with the canonical distribution at a fixed 
value of X] the corresponding free energy is given by the expression 

F{X) = -(3-HnZ{X). (13) 

Combining (|TOl) . (|TT1) . and (fT3ll we get, for a conjugate pair of trajectories 7^ and 7-^, 

where AF = F{B) — -F(A). We now use (HD, ([S]), and ([H]) to rewrite the quantity inside 
the last exponent above: 

H{z^; B) - H{z^- A) = H{z^*- B) - Hiz^; A) 

= H{z^- B) - H{z^- A) 

= W, (15) 

where W represents the work (as defined by ([8])) performed on the system during the 
realization 7^. Thus we finally obtain 



Pf[i 



F^ 



,,iw-.n_ (16) 

This result was first derived by Crooks [8], who modeled the evolution of the system 
as a stochastic discrete time Markov process, rather than using Hamiltonian dynamics. 
One consequence of (fT6l) is the Crooks Fluctuation Theorem for the ratio of forward 
and reverse work distributions [9]: 

Pf{W) _ ^I3^w-AF) 



mw-/\F) 

Pr{-W) ^ ■ ^'^) 

Often (ITBl) and ffT7|) are expressed in terms of the dissipated work, = W — AF ^5] . 

2.3. Bochkov-Kuzovlev Fluctuation Relation 

In deriving f|T6|) . we imagined that initial conditions were sampled from the canonical 
distributions associated with the Hamiltonians H{z]A) and H{z;B), for the forward 
and reverse processes, respectively. Now let us consider a different situation, in which 
" for both processes - initial conditions are sampled from the canonical distributions 
corresponding to the bare Hamiltonian Hq{z). We suppose, however, that the protocols 
X^ and X^ are the same as above, namely from A io B and from i? to A, where 
these values remain arbitrary. If we use Pp[7'^] and -Pr[7^] to denote the corresponding 
distributions of trajectories, then we have 

^F[7l = ^exp[-/3/7o(^o^)], (18) 

i^M = ^exp[-/3/Jo(^o'^)] (19) 
in place of ffTOl) . ffTTj) . Taking the ratio then gives us 
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By analogy with (IT^ . we have 

Ho{z^) - Ho{z^) = Ho{z^) - Ho{z^) = Wo, (21) 

where Wq is the work (as defined by ([7])) performed during the reahzation 7^. Thus we 
finally arrive at 

which corresponds to equation 2.11 of [Ij. 

Related to (!22|) is a fiuctuation theorem for the work distributions (cf. (fTTI) ). 

P°f(M'o) _j,Wo 



where Pp/jiiWo) are the work distributions associated with the trajectory distributions 
in f|T8l) and f|T9|) . Note that the quantity Wq that appears in fl22|) and fl23l) is in general 
different from the dissipated work Wd- 



As highlighted by these derivations, a principal difference between ( fT6l) and ( l22l) 
(equivalently, and ([2])) is the choice of the equilibrium distributions from which 
initial conditions are drawn. If we sample from canonical distributions corresponding 
to the Hamiltonians H{A) and H{B), then we obtain f|T6|) . while if we sample initial 
conditions from the bare canonical distribution, corresponding to Hq, then we arrive at 
(122]). 

It is useful to spell out the physical interpretation associated with this distinction. 
In Crooks's formulation (JT6ll . we imagine that prior to t = we allow the system to 
come to equilibrium with a thermal reservoir at inverse temperature /3, while holding the 
external force fixed at either A or B. We then let the system evolve as the external force 
is manipulated according to the forward or reverse protocol. To obtain Bochkov and 
Kuzovlev's result (!22|) . we instead imagine that the external force is absent (X = 0) as 
the system equilibrates with the reservoir. The system thus settles into an equilibrium 
state corresponding to the bare Hamiltonian Hq. Just prior to t = we suddenly turn 
the force on to A or 5, subsequently switching it continuously according to the forward 
or reverse protocol. 

The two results involve not only different distributions of initial conditions, but also 
different definitions of work, Wq and W. The right sides of (JTj) and ([H]) suggest a physical 
interpretation of this distinction. Namely, in one case we define the internal energy of 
the system to be given by the value of the bare Hamiltonian Hq; then Wq is viewed as 
the work performed by the application of an external force X that affects the evolution 
of the system in a fixed energy landscape Hq. In the other case, the system's internal 
energy is taken to be the full Hamiltonian H; we then view W as the work associated 
with externally controlled manipulations of the energy landscape H = Hq — Xa. See 
[TT] for a more detailed discussion, and illustration, of this distinction. 

In the special case of cyclic processes, defined by the condition Xq = X^, we have 
AF = and W = Wq [TT]. In this situation (fT6l) and (l22ll are equivalent. If we further 
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restrict ourselves to time symmetric protocols (Xj = X^_t), then there is no distinction 
between the forward and reverse processes. In this case, f|T6|) and fl22l) can be viewed as 
special cases of a more general result derived by Evans and Searles (see equation (4.15) 
of [21], as well as the recent review by Sevick et al [25j). 



3. Extended Relation 



In section [2] we used Hamiltonian djTiamics to model the evolution of the system. By 
contrast, as noted above. Crooks modeled this evolution as a discrete-time, stochastic 
Markov process, to represent a system in contact with a thermal reservoir [8|. If we were 
to carry out an analysis similar to that of section [2], but in a Markovian rather than a 
Hamiltonian framework, then we would arrive at the same conclusion, namely that both 
([T]) and (l2l) are valid, the difference between them relating to the choice of canonical 
ensemble from which to sample initial conditions. Rather than showing this explicitly, 
in this section we derive a more general result, (l36l) below, that contains both ([1]) and 
([2]) as special cases. We will present both Hamiltonian and Markovian derivations of 
this result. 

Consider a Hamiltonian of the form 

H{z; Xt, Xt) = Ho{z; Xt) - Xta{z), (24) 

where now Hq itself is parametrized by a variable A. Thus the full Hamiltonian H is 
made time-dependent by manipulating both the force X, and a work parameter, X. From 
a mathematical perspective, this division of the parameters into an external force and 
an internal work parameter is somewhat arbitrary, particularly if A appears linearly in 
Hq (as in ( H9l) below). While in specific situations the distinction may be motivated by 
physical considerations or by the questions being addressed, our aim here is to develop 
a general formalism that encompasses both ([T]) and (I2l). For this purpose, the form of 
H given in flM|) is useful. 

For the bare Hamiltonian Hq, we have a parameter-dependent canonical distribution 

Pl\z;X) = ^^e-^^^'--''\ (25) 

for which the partition function and free energy are given by 

Z(A) = y dze-^^°("'^) and F{X) = -p'HnZ{X). (26) 

Throughout this section we will take the internal energy of the system to be defined 
by the value of the bare Hamiltonian Hq, and we will consider processes during which 
the system evolves in time as both A and X are switched according to pre-determined 
protocols. Thus we picture a system that evolves in a time- dependent energy landscape 
Hq, while also coupled to a time-dependent external force X. For a given realization of 
such a process, we take the work performed on the system to be given by 

riH 

W= / A— dt+ / Xadt, (27) 

./n (jX In 
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or by the discrete-time analogue of this expression, (1391) below. This is a hybrid definition 
whose two terms are similar to the expressions for W and Wq introduced earlier. As 
we will see, for processes during which the system is thermodynamically isolated the 
value of W is equal to the net change in Hq, (133|1 : whereas if the system remains in 
contact with a reservoir as A and X are varied, then the net change in Hq is equal 
to the sum of W and a quantity Q that represents the heat absorbed by the system, 
f l42|) . In either case, fl271) provides a definition of work that is faithful to the first law of 
thermodynamics, and consistent with our definition of internal energy, Hq. 

As before, we will compare two processes (F, R), characterized by conjugate 
protocols for the externally controlled parameters: 

^/^ = ^f-. , Af = Af_,. (28) 

For the forward process, we let A and B denote the initial and final values of X, as 
before, and we let and 1 denote the initial and final values of A. For the reverse 
process, then, X is switched from B to A, and A from 1 to 0. We will assume that 
initial conditions for the forward and reverse processes are sampled from p^{z] A = 0) 
and p'^{z; A = 1), respectively (see fl25|l ). The physical interpretation of this assumption 
is similar to that discussed in the context of the Bochkov-Kuzovlev results, near the end 
of section [21 Namely, we imagine that the external force is absent (X = 0) as the system 
equilibrates with the reservoir prior to time t = 0. The work parameter is held fixed 
at either A = (for the forward process) or A = 1 (for the reverse process) during this 
equilibration stage. Then, immediately prior to t = the external force is turned on to 
either A or B, and subsequently both A and X are varied according to the appropriate 
protocol. 



3.1. Hamiltonian dynamics 

With the setup just described, we now suppose that the microscopic history of the system 
during the time interval < t < r (for either process) is described by a Hamiltonian 
trajectory Zf. Such evolution has the property that the total rate of change of the 
Hamiltonian is equal to the partial derivative of the Hamiltonian function with respect 
to time (see e.g. equations 8-35 of [22]). Thus 
d f)M f)M 

-/f(.,;A„A,)^A,- + A-,- (29) 

= \^-Xa (30) 
oX 

f)J-f d 

= X,^ + Xa--iXa). (31) 

Adding (d/dt)(Xa) to both sides, we get 
d flM 

_ff„ = A._^ + .«, (32) 

The first two terms on the right are the integrands appearing in our definition of W 
([27]). Thus, integrating both sides of (1521) with respect to time, we see that W is the 
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net change in the value of the bare Hamiltonian: 

W = Hoiz,;Xr)-Hoizo;Xo). (33) 

This identity will now be put to use much as and ([H]) were in section O 

Adopting notation similar to that of section [21 we let -Pf[7^] and -Pr,[7'^] denote the 
probability distributions of realizations of the forward and reverse processes; as before 
we use the Liouville measure of initial conditions to define a measure on the space of 
possible trajectories, for both processes. We then have, by analogy with ( fTOl) and ( fTTl) . 

Ppb^] = pT(4; 0) = ^ exp[-l3Ho{z^; 0)] (34) 

Pnb^] = f^{z^- 1) = ^ eM-f^Mz^: 1)], (35) 

since initial conditions are sampled from the canonical ensembles associated with Hq. 
Taking the ratio of these two distributions for a conjugate pair of trajectories, and 
making use of fl26|) and fl33|) . we obtain 



Pf[i 



^m-AF)^ (36) 



where AF = F{1) - F{0). 



3.2. Markovian Stochastic Dynamics 

We now derive the same result, 0361) . for a system that remains in weak contact with a 
thermal reservoir during the interval < t < t, and whose evolution in this situation is 
modeled as a discrete-time Markov process. Specifically, we imagine that the microscopic 
history of the system is described by a sequence zq, Zi, ■ ■ ■ zjy, representing phase space 
points visited at times to^^i) ■ ■ 'tN, with tn = ut/N. With each time increment, both 
the phase space point z and the parameters vr = (A, X) are updated with a combination 
of a stochastic step and a switching step, as follows. 

Stochastic step. During a stochastic step, starting from, say, Zn at time tn, the 
system makes a random jump to a new point in phase according to the 

transition probability P{zn Zn+i;X,X), where (A,X) are the fixed values of the 
external parameters during this step. This transition probability is assumed to satisfy 
detailed balance (a requirement that may be reduced, at the expense of using a modified 
dynamics for the reverse process [23]), expressed as 

Pizn^Zn+i;X,X) ^ e-^^(-"+^^^-^) 

P{zn^z^+,;X,X) e-/5^(^";^.^) • ^'"^ 
The direction of the arrow indicates the sense of the transition, e.g. from Zn+i to Zn in 
the denominator. Since each stochastic step, z„ —* -^n+i, is meant to represent effects 
of the interaction with the thermal reservoir, any change in the system's energy during 
such a step will be interpreted as heat exchanged with the reservoir; see fl55]) below. 

Detailed balance (1371) represents an assumption about the dynamics of the system. 
Among other things it implies that the equilibrium distribution is conserved when the 
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parameters A and X are held fixed. Since the equilibrium distribution is determined 
by all the physical forces acting on the system, both those that we choose to view as 
internal and those that we treat as external, detailed balance is expressed in terms of 
the full Hamiltonian, H, rather than Hq. 

Switching step. During a switching step the external parameters are updated, e.g. 
from Tin = (A„,X„) to Hn+i = (A„+i,X„+i), according to a given protocol, with the 
phase space point held fixed. 

The evolution of the system from to = to t at = r consists of an alternating 
sequence of stochastic and switching steps. For the forward process we assume that the 
switching step comes first, giving us a progression of the form {z^, n^) {z^, tt^+i) —>■ 
{z^+i,TT^j^_i) at each time step. For the reverse process, the stochastic steps leads: 

We now define the heat absorbed by the system from the reservoir, for a given 
realization of the forward process, as a sum of energy changes during the stochastic 
steps: 

N-1 N-1 

Q = ^ SQn = ^ [H{Zn+i; TTn+l) " H^Zn, TTn+i)] . (38) 
n=0 n=0 

For the same realization, the work performed on the system is defined as the discrete- 
time version of ( 1271) : 

Af-l 

W= J2^^n 
n=0 

7V-1 Af-1 

= ^ [Ho{zn; A„+i) - Ho{zn; A„)] + ^ [(a„+i - a„)X„+i], (39) 

n=0 n=0 

introducing the convenient shorthand a„ = a{zn)- On the right side of (l39l) . the first sum 
represents energy changes due to the switching of A, while the second sum is associated 
with the force X acting on the coordinate a. From these definitions, the identity 

6Hn = H{Zn+l] TXn+l) " H{Zn] TTn) 

= 5Qn + SWn + anXn - a„+iX„+i (40) 

is easily verified, by inspection. Combining fl38l) . fl39l) and fj40l) now gives us the following 
expressions for the net changes in H and Hq, from t = and t = r: 

H{zn; tin) - H{zo; no) = Q + W + XqUq - XnUn (41) 
Hq{zn; \n) - Hoizo; Xo) = Q + W. (42) 

The latter is analogous to fl33|) . but now includes heat exchange with the reservoir. 

With these preliminaries out of the way, we analyze the probability distributions 
of trajectories, for the forward and reverse processes. A trajectory is specified as a 
sequence of + 1 points in phase space, e.g. 7 = {zo,Zi, ■ ■ ■ z^)- A natural measure 
on the space of such trajectories is given by d7 = dzodzi ■■■dzN- The probability 
of obtaining a particular trajectory, P['j], is equal to the probability of samphng its 
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initial conditions zq from the appropriate equilibrium distribution, multiplied by the 
conditional probability of generating the subsequent sequence of steps, -P[7|zo]- Thus 
for the forward process we have 

^^[7^ = Ao^ = 0) ■ PAl^^] (43) 

PAi^\4] = P{'4 - 4; <) ■ --Pi^-i - 4; 4) (44) 

For the reverse process we have analogous equations, except that = 1, and the 
transition probabilities are of the form P{z^ ^n+i] '^n)y reflecting the convention that 
the stochastic step precedes the switching step, during each time interval of the reverse 
process. In writing the conditional probability as the product of transition rates in (l44ll . 
we have made use of the assumed Markov nature of the dynamics. 

For a conjugate pair of trajectories 7^ and 7*^, the ratio of conditional probabilities 

is: 

Prb^lz^] ^ P(^o^->zf;vrf)---P(z^„,->z^;7r^) 

pRb^K] P(^o^-^f;vrf)---P(z^_,-z^;7r^)' ^ ^ 

where we have used the relations = z^_^ and = 7r^_„ to rewrite Pjib^lz^] in 
terms of z^'s and vr^'s. Following Crooks [8], we now use detailed balance to express 
the right side of fHSj) as a product of factors exp{—(3SQn) (see fl37j) . fl38|) ). which gives 
us 

1^1^ = e-^<^. (46) 
PrIi'^K] 

To get the ratio of unconditional probabilities Pp[7'^]/Pr[7'^], we multiply the numerator 
and denominator of both sides by the appropriate initial distributions, as per (j^3l) : 



using Zq = zff. With (j42j) this becomes 



Pij[7^] 

as in the Hamiltonian case (!36l) 



3.3. Discussion 

We now argue that ([1]) and ([2]) can be viewed as special cases of fl36l) . To establish 
this, let us suppose that the parameter A appears linearly in Hq, so that the full, time- 
dependent Hamiltonian takes the form 

H{z; Xt, Xt) = Hoiz) - Xt viz) - Xt a{z). (49) 

It is natural to view both A and X as forces, coupled to the coordinates rj and a, 
respectively. However, in keeping with the convention adopted throughout section [3l we 
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treat the term —Xrj as a contribution to the internal energy of the system, Hq = TYq — A?7, 
and X as an external force. The work, fl27|) . is now given by 

(50) 

With these 



W = - \ridt+ / Xadt, 
Jo Jo 

or by the corresponding discrete-time analogue obtained from fl39l) . 
definitions ( l36i) becomes, exphcitly. 



PrIi 



R] 



exp Ar/dt + y Xddt-AFj . (51) 

As before, initial conditions are sampled from the canonical ensemble p'^ oc e~^^°, 
associated with the initial value of A (A^ = and A^ = 1). 

If we now consider the special case in which A = during the entire process (hence 
oc e~^^° and AF = 0), then the right side of (15T1) becomes exp (/? / Xce dt) , and we 
recover ([2]). If we instead imagine that X = during the entire process, then the right 



(3 i — J Xrj dt — AF ) . We then recover ([T]), with A and ri{z) 



side of ( 1511) becomes exp , , j , 
playing the roles assigned to X and a{z) in section [221 

To summarize, if we view ([1]) as a result that applies when a system evolves in a 
time-dependent energy landscape, and ([2]) as a result that applies when a system evolves 
under the influence of a time-dependent external force (as per the discussion at the end 
of section [2]), then (!36|) is a generalization that applies when both factors are present. 



4. Conclusion 



Our central aim has been to clarify the relationship between the fluctuation relation 
obtained by Crooks, ([1]), and that due to Bochkov and Kuzovlev, ([2]). These two results 
share a similar structure, but are not equivalent. The distinction between them is traced 
to the different definitions of work {W, Wq) used by the authors. One consequence of this 
difference is that ([1]) involves (and can be used to estimate) a free energy difference, AF, 
while (I2l) does not. Moreover, the two results are associated with different definitions 
of the internal energy of the system, and rely on different assumptions regarding the 
initial equilibrium preparation of the system. Finally, we have obtained a generalized 
fluctuation relation fl36l) that unites ([1]) and ([2]). 
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